ITEP-TH-34/98 
|hep-th/9806T50| 
June, 1998 



Interquark Potential in Schrodinger Representation 



K. Zarembo 



Department of Physics and Astronomy, 
University of British Columbia, 
6224 Agricultural Road, Vancouver, B.C. Canada V6T 1Z1 

and 

Institute of Theoretical and Experimental Physics, 
B. Cheremushkinskaya 25, 117259 Moscow, Russia 

E-mail: zarembo@itep . ru/@theory .physics .ubc . ca 



Abstract 

Static charges are introduced in Yang-Mills theory via coupling to heavy fermions. 
The states containing static color charges are constructed using integration over gauge 
transformations. A functional representation for interquark potential is obtained. This 
representation provides a simple criterion for confinement. 



1 Introduction 



The gluon vacuum is expected to determine many properties of the low-energy QCD. Ac- 
cording to the most straightforward definition, the vacuum is described by a ground state 
wave function - the lowest energy solution of the stationary Schrodinger equation. The exact 
ground state in gluodynamics must be very complicated, but one may hope that a reasonably 
good approximation can be obtained by means of variational or of some other approximate 
method if a proper parametrization of the wave functional is used. Such parametrization 
should respect the symmetries of the Yang-Mills theory, first of all, the gauge invariance. 
There is a simple way to achieve the invariance of a wave functional - to enforce it by 
integrating over all gauge transformations [2|, [| f| : 

V[A}= J[DU}e- s[AU \ (1.1) 

A? = U^A + ^d)jU. (1.2) 

Alternative approaches are based on the use of gauge-invariant variables from the very be- 
ginning || or on fixing of residual gauge symmetry ||. 

The representation of a projector on the physical states in a form of the integral over 
gauge transformations is well known in Yang-Mills theory at finite temperature J/J. In the 
present context this method appears to be very convenient for variational calculations, since 
it allows to make the simplest, Gaussian variational ansatz JI], [], |j] exactly gauge invariant. 

The averaging over the gauge group introduces additional variables, the parameters of 
gauge transformations, which more or less correspond to longitudinal components of the 
gluon fields. Introduction of such variables might seem somewhat artificial, since longitu- 
dinal degrees of freedom do not propogate and, thus, are unphysical in pure gauge theory. 
This is true, however, until external color charges are considered. Electric fields of a static 
charge are longitudinal, so longitudinal components of gauge fields are responsible for color 
electric forces, which, in their turn, are responsible for confinement. Because averaging over 
gauge transformations explicitly separates longitudinal degrees of freedom, it provides an 
appropriate framework for consideration of static color charges and of their interaction ||. 
In this paper the averaging over gauges is applied to Yang-Mills theory coupled to heavy 
fermions, which allows to introduce static charges in a more systematic way. 



2 Coupling to heavy fermions 

The ground state in gluodynamics satisfies Schrodinger equation 

H YM V = E vac V, (2.1) 
where Hym is Yang- Mills Hamiltonian (in A = gauge): 



H\m= !<?■>■ (\EfEf + \f$F£) . (2.2) 
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When the wave function is a functional of the gauge potentials, the electric fields act as 
variational derivatives: 

E <—'W < 2 - 3 > 

We consider SU(N) gauge group with generators T A normalized by tr T A T B = —5 AB /2 and 
sometimes use matrix notations for gauge potentials: A4 = A A T A . 

The Hamiltonian (2.2) commutes with operators DiE A , and the ground state wave func- 
tional is subject to the Gauss' law constraint: 

DiEf ^ = 0. (2.4) 

This constraint represents in the infinitesimal form the gauge invariance of the ground stateQ: 

y[A n }=y[A]. (2.5) 

The integral (|1 . 1|) can be considered as a general solution of the Gauss' law constraint. Any 
gauge- invariant state can be represented in this form, the exact ground state in gluodynamics, 
in particular. Of course, the representation ( |1 . 1| ) is aimed to the use of approximate methods, 
when some simple form of the action S [A] is employed. In the variational approach of Ref. |[| 
the action was taken quaratic and diagonal in spatial indices. This variational ansatz was 
generalized to include longitudinal part of the quadratic form [|J and tree-level condensates 
H . However, since the consideration below is purely kinamatical, we do not restrict ourselves 
to any particular form of the bare wave functional exp(— 5*[A]). 

In order to introduce color charges, we couple the gauge fields to heavy fermions, so 
heavy that it is possible to neglect their kinetic term in the Hamiltonian: 

H = J d 3 x Q E A E A + 1 F A F A + . (2.6) 

The fermion operators obey anticommutation relations 

{r a (x),4 b (y)} = 5 aP 5 a b 8(x-y), (2.7) 

where a,b and a,/3 are color and spinor indices, respectively. The fermion fields transform 
in the fundamental representation of SU(N): 

ip u = tfip, ^ U = ^U. (2.8) 

Physical states are gauge invariant and satisfy the Gauss' law constraint: 

(D t E A + % ^T A i) # phys = 0. (2.9) 

Because of the absence of the kinetic term for matter fields, fermion operators obey 
simple commutation relations with the Hamiltonian. In the standart representation of the 
Dirac matrices, when 7 is diagonal: 

7°=( \ _°A (2.10) 



"For simplicity, we disregard peculiarities related to topologically nontrivial gauge transformations, see 
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the commutation relations take the form: 

[H,i/> a (x)] = ±Mi{> a (x), ZZi,$ ( 211 ) 

H,^l(x)} = ±M£(x), I +| • (2-12) 

Therefore, ^3,4(2;), ip[ 2 (a;) and ^i )2 (a;),"03 4(^) are creation and annihilation operators, re- 
spectively. The fermion vacuum is defined by equations 

Vi, 2 (a;)|0) = 0, 

^J )4 (ar)|0) = 0. (2.13) 

The fermion vacuum is gauge invariant by itself, and the ground state of the Hamiltonian 
( |2.6| ) is described by the wave function 

vf vac = *[A]|0) = J [DU] e- s ^ |0>, (2.14) 

where *&[A] is the vacuum wave functional in gluodynamics - the lowest-energy solution of 
the Schrodinger equation ( [2.1| ). 

Excited fermion states are obtained by acting of creation operators on the vacuum. These 
states are no more gauge invariant, and it is necessary to average them over gauge transfro- 
mation in order to project on the physical subspace. The physical states are, consequently, 
constructed with the help of dressed operators 

€ a \v) = u^'SvWM, 

4%(x) = 4 b (x)U b b ,(x), (2.15) 

where a = 3, 4 and (3 = 1, 2 (below spinor indices are omitted for brevity). The most general 
state of this type, 

K> '•<■•//, »,) = f[DU] e- s ^4 a F(x 1 )...4 a F(x n ) 

i> uv i(yi)...i> u HysM, (2-16) 

describes n static antiquarks and s quarks located at points Xi, yi. Dressed operators ( [2.15| ) 
are invariant under local gauge transformations because gauge rotations of the gluon and of 
the quark fields; A± — > Af, ip — ► ip , ijj* — > ip* ] are compensated by a change of variables 
U a a , -> Q U b U b a , in the integral over the gauge group. These transfromations do not act 
on indices with prime. To be more presize, this statement is true for transformations with 
lim Q(x) = 1. If lim Q(x) = Qoo 4 1, the change of variables U — > Q*U violates boundary 

x— KX x— >oo 

conditions in the path integral over U and should be accomponied by a global transformation 
U — > UQoo. This means that indices with prime are global color ones. 

The state ([2.16|) can be factorized in a product of the gluon and the fermion wave func- 
tions. For example, quark-antiquark pair in a color-singlet state is described by the wave 
function 

<M*,y) = / [DU] e- s ^4F(x)i; Ua '(y)\0) = *\[A; x, y] 4 a (x)^ b (y)\0), (2.17) 
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where 

* a b [A;x,y) = J [DU] e~ s ^ U a a ,(x)tf a ' b (y). (2.18) 

The gluon wave function (|2.18| ) transforms in the representation N <8> N of the gauge group: 

* a b [A n ;x,y] = rt a c (x)n d b (y)* c d [A;x,y] (2.19) 

and satisfies Schrodinger equation for pure Yang-Mills theory. 

Again, any state obeying transformation law ( |2.19| ) can be represented in the form of 
an integral ( |2.18| ); moreover, the same action S[A] can be used both in the vacuum and 



in the charged sectors. The integration over gauge transformations projects the bare wave 
functional exp(— S[A\) on subspaces of states with definite transformation properties, gauge 
invariant states in the former case and the ones belonging to N £g> iV representation in the 
latter. These subspaces are orthogonal in the full Hilbert space, and it is always possible 
to find a state which has given projections on each of them. Hence, the action S[A] can 
be chosen in such a way that the projections of the wave functional exp(— S[A]) on the 
sectors defined by transfromation laws ( |2.5| ) and ( 2.19p are the lowest-energy eigenstates of 



the Yang-Mills Hamiltonian in each sector. So, the vacuum, the quark-antiquark and, in 
fact, any state of type ( |2.16 ) can be obtained by averaging of the same bare wave functional 



over gauge transformations. The above arguments concern exact wave functionals, but the 
statement that one can use the same action S[A] in different sectors is also true within a 
variational approach. Suppose that the form of S[A] is fixed up to some parameters, then 
variational estimates for these parameters will be the same in the vacuum and in the charged 
sectors because energies of the ground states in these sectors, being proportional to the 
volume, differ by an amount which is finite in the infinite volume limit and, thus, does alter 
variational equations. 

Another example of the state ( 2.16|) is the one with baryonic quantum numbers: 



<Ml/i,...,Jfor) = J[DU]e- s ^ U h a[ ^J u <(y 1 )...r a 'HyN)\0) 

= [4 yi,...,y N ] VKvi) ■ ■ ■ rHVNM, (2.20) 

^ai...a N [A; 2/i, ...,Vn] = J[DU] e"^ e^U^Jy,) . . . U U '\ n (y N ). (2.21) 

The state (|2.18|) has simple physical meaning, it describes color fields of static quark- 
antiquark pair in the gluon vacuum. In the Abelian theory, Coulomb field of a static charge 
is created by an operator e %eV ( x ) where 

V(x) = J d 3 y - y)diAi(y), (2.22) 

so dressed electron operators have the form ||: 

^:\y) = e- ieViy) My), 

^* ] \x) = e ieV( - x ^l(x). (2.23) 

It is not difficult to reproduce this construction from the representation based on averaging 
over gauge transformations. In Abelian theory U = e ie(p , Af = Ai + di^p and the action S[A] 
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should be quadratic; for simplicity we choose the diagonal quadratic form. The integration 
over gauge transformations in (|2.18|) is Gaussian and can be done explicitly: 



*[A;x,y] 



[dip] exp 



--( d 3 xd 3 w {Ai + dup) (z) K{z - w) (Ai + di<p) (w) 



+ ie<p(x) — ie<p(y) 
C e ieV{x) e- ieV{y) exp 



l - J d 3 zd 3 wA^{z)K{z - w)Al(w)\ , (2.24) 



where C is a field-independent constant and Aj~(z) denotes transversal part of the gauge 
potentials: 



A 




(2.25) 



As follows from conformal symmetry, the coefficient function K is equal to \p\ in the momen- 
tum space. Then the last factor in ( |2.24| ) is nothing but the vacuum wave functional for free 
electro-magnetic field. The first two factors reproduce dressing operators which correspond 
to a pair of static charges of opposite sign located at x and y. More generally, the integration 
over the Abelian gauge group in ( |2.16|) , due to its Gaussian nature, leads to the replacement 
of ip u by ip(*\ so dressing ( |2.15| ), ( |2.16| ) is equivalent to ( |2.23| ) in the Abelian theory. 

Dressed operators ( |2.23| ) possess a number of important properties. They are gauge- 



invariant and create eigenstates of the free electro-magnetic Hamiltonian. In QED they also 
play an important role, because dressed electron and positron operators allow to avoid IR 
divergencies related to emission of soft photons [II], [II]]. Different non- Abelian generalaza- 
tions of the operators ( |2.23| ) were proposed |I| [L2|]. In contrast with these constructions, 
dressing (|2.16| ) in non-Abelian case does not have factorized form, the integral over gauge 
transformations is no more Gaussian, and the operator creating static gauge fields of one 
particular charge cannot be defined; that is, the form of the wave functional describing this 
charge depends on the presence of other charges. 



3 Interquark potential 



The energy of the state (|2.17|) after obvious subtractions determines qq interaction potential: 

2M + E^ c + V{x-y). (3.1) 





H\V M ) 







Matrix elements of gauge- invariant states, like those entering eq. (|3.1|) , are proportional to 
the volume of the gauge group. The representation ( |1.1| ) allows to get rid of this infinite 
factor easily |IJ. For example, the norm of the vacuum state is 

(#|#)=/ [DU][DU'][dA] e s[A u ]-s[A^] = congt J e -S[A«]-S[A}_ (3 _ 2) 
The norm of the charged state ( |2.17| ) is 



(y M \V M ) = const J [DU][dA] trU(x)U^(y)e- s[AU] - s[A] . (3.3) 
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Matrix elements of the Hamiltonian correspond to insertion of the operator R[A], 



1 5 2 S[A] 1 5S[A] 5S[A] 1 



if variational derivatives in the Hamiltonian act on the right. If they act on the left, R[A] 
is replaced by -RLA^ 7 ]. It is convenient to introduce a source for the symmetric combination 
(R[A u }+R[A})/2: 

~[DU][dA] e-^l-W^A^BM (3.5) 

then matrix elements of the Hamiltonian are obtained by differentiation with respect to A. 
For example, the vacuum energy is given by 



E vac = — — In Z 

OA 



(3.6) 

A=0 



This quantity is proportional to the volume and contains UV divergent contribution from 
zero-point oscillations. Regularized energy density can be related to gluon condensate [|13| : 



16a s \- 

where /3(oe s ) is a /3-function in Yang-Mills theory. 
For qq potential we obtain: 

V(x-y) = -^\n(trU(x)U\y) 



0) Vol + UV divergent terms, (3.7) 



(3.8) 

A=0 

where averaging is defined by the partition function (|3.5|) . It is worth explaining how con- 
finement can arise in such representation. The confining potential grows with distance; at 
the same time, the correlator in ( j3.8|) should decrease at infinity. But its logarithm can 
increase thus leading to the confining potential. It happens when the integration over the 
gauge group produces a mass gap: 

(tiU{x)U\y)) oc e~ mr , (3.9) 

where r = \x — y\ — > oo. Then interquark potential grows linearly with distance: 

V(r) = ar + ..., (3.10) 

the string tension being equal to the derivative of the mass: 

dm 



(3.11) 



A=0 



Hence, the confinement arises due to generation of a mass gap in the averaging over gauge 
transformations. More presizely, the confining potential is generated if there is a nonzero 
linear response of the mass gap to switching of the Hamiltonian. 
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Analogously, the energy of quarks in the baryonic state fl2.2U|) , 



V B (yi, ■ ■ -,Vn) 







(*B 


#b) 



NM -E v 



(3.12) 



can be represented as the derivative of certain correlation function in the statistical system 
defined by the partition function (|3.5|): 



VbGfr,...,!/*) = ^ln(e ai ... aiV ^- fe -f/ t \( Z/l )...^ ta V(^) 



(3.13) 



A=0 



4 Discussion 



Examples of calculations of interquark potential starting from representation based on av- 
eraging over gauge transformations are given in [|2|, [|. These calculations use particular 
assumptions about the structure of the ground state wave functional in gluodynamics and 
are by no means derivations from first principles. On the level of kinematical consideration, 
one can state that the potential is confining if the integration over gauge transformations 
generates a mass gap with nonzero linear response on switching of the Hamiltonian. The 
calculations [fj, |4[] show that generation of a mass gap is rather natural in non-Abelian the- 
ory, but quantitative results can be obtained only if the vacuum wave functional is known in 
some approximation. Perhaps, the simplest Gaussian ansatz averaged over gauge transfor- 
mations m, H|, [| can provide a reasonable variational approximation. The main problem is 
correct renormalization. The variational method always overestimates the vacuum energy, if 
non-renormalized divergencies remain, the overestimate is infinite, and variational approxi- 
mation loses sense [13 . It was argued M that Gaussian ansatz can be made compatible with 
asymptotic freedom; the renormalization properties of Yang-Mills theory in the Schrodinger 



representation were also discussed in [15] from different points of view 
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